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Abstract 



We present a method how to detect the W ± H° — > l^vbb in the high luminosity LHC environment 
with the CMS detector. This study is performed with fast detector response simulation including high 
luminosity event pile up. The main aspects of reconstruction are pile up jet rejection, identification of 
fe-jets and improvement of Higgs mass resolution. 

The detection potential in the SM for m H a < 130 GeV/c 2 and in the MSSM is only encouraging for 
high integrated luminosity. Nevertheless it is possible to extract important Higgs parameters which 
are useful to elucidate the nature of the Higgs sector. In combination with other channels, this channel 
provides valuable information on Higgs boson couplings. 
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1 Introduction 



The observation of at least one Higgs boson is an important proof of the Higgs mechanism |[TJ which is introduced 
to explain the masses of elementary particles. Beside the discovery of a Higgs boson, it is also important to study 
the Higgs boson couplings. Higgs bosons lighter than 130 GeV/c 2 decay mainly to a bb pair [ph. In this note we 
describe a method to observe a Higgs boson in the associated production channel W^H — > l^ubb, shown in 
Figure [l], with leptonically decaying W ± . Among other production channels with H° — > bb decay [||, the W ± H° 
channel turns out to have a low signal to background ratio and therefore requires a large integrated luminosity. The 




Figure 1: W^H — > V^vbb signal event at LO. 



signal and backgrounds are simulated in the high luminosity regime with each event superimposed on an average 
of 17.3 Poisson distributed pile up events generated with PYTHIA [§] (MSTP(2) = 2 , MSTP(33) = 3 , MSTP(81) 
= 1 , MSTP(82) = 4 and PARP(82) = 3.). The relevant signal and background cross sections at the LHC (y/s pp — 
14 TeV) and particle masses used in the simulation are listed in Table |lj 



LO cross sections 


masses 


<rw±H% M x BR H « M ^ b - b = 2.51 - 0.41 pb 
a w ± z<> = 27 pb 
<Jw±jj = 30 nb 
on = 570 pb 
<j tb - = 320 pb 


m H o = 90- 135 GeV/c 2 
m z o = 91.187 GeV/c 2 
m w ± = 80.41 GeV/c 2 
m t = 173.8 GeV/c 2 
m b = 4.3 GeV/c 2 



Table 1: PYTH1 A cross sections of signal and background relevant for the W H° — > l^vbb channel and calculated 
with parton density function CTEQ41 jj^. The branching ratio of the decay to electrons ormuons (= 22%) is 
not included in the cross sections of this table. The particle masses are from /HJ. The generation of single top (tb ) 
events can be done in PYTHIA with processes number 2 and 83 by forcing v§ > 180 GeV jj^. 

After the event generation, the detector simulation is performed: FATSIM [^J is used for the simulation of tracker 
response and track reconstruction. Track momentum smearing, impact parameter smearing, impact parameter tails 
and track reconstruction efficiency are taken into account, as well as geometrical acceptances. CMS JET is used 
for the simulation of calorimeter response, jet reconstruction, missing transverse energy calculation and lepton 
smearing. The reconstruction efficiency for leptons (in this study only electrons and muons) is assumed to be 90%. 
All parametrisations have been obtained from GEANT based detailed simulations. 

2 Reconstruction 

A typical signal event as shown in Figure [j] is expected to give a final state consisting of one isolated lepton, 
missing transverse energy E™ lss and two 6-jets. Pile up events, underlying event and gluon radiation are sources 
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of additional jets which complicate the reconstruction and selection. An initial study JIo| ] of this channel showed 
already the importance of two reconstruction features, 6-tagging performance and jet reconstruction: 

A clean and efficient identification of 6-jets is important to reduce most of the background with less than two 
genuine &-jets. For a realistic (and reasonably fast) 6-tagging simulation, the understanding of impact parameters, 
impact parameter errors and track reconstruction efficiency are crucial. In Figure [| the parametrisation FATSIM 
(fast tracker simulation) is compared with the detailed simulation: the ^-tagging performance of the CMS tracker 
(phase 1) is simulated for soft jets coming from the Z° decay using mainly the significance of the signed transverse 
impact parameter a(ip) cuts on two tracks per jet. There is a good agreement between the two simulations. In 
this example, the 6-tagging efficiency for a(ip) > 2 is 50% and the mistagging probability is less than 1%. The 
6-tagging efficiency for jets from the Higgs decay is higher (« 60%), because the jets coming from the heavier 
Higgs boson are more energetically. Identification of leptons inside jets and reconstruction of secondary vertices 
can improve the 6-tagging performance further. 
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Figure 2: b-tagging performances of Z° — > jj events by FATSIM (bright) and CMSIM (dark). The tagging rates of 
b- and u-jets were determined with the following algorithm: 2 tracks in a cone of 0.4 around the jet (Et > 20 GeV, 
\t]\ < 2.4) axis with pr > 0.9 GeV/c, min. 6 hits including 2 pixel hits (in the case of CMSIM) , ip < 2 mm and 
the signed transverse impact parameter significance cr{ip) > 0.0 ... 4.5 (various cuts). 



A good Higgs mass resolution helps to enhance the signal event invariant m(j,j) mass peak allowing a better 
visibility of the signal and a more precise determination of the Higgs mass. A good mass resolution can be 
obtained, when the energy and direction of each reconstructed jet agree as closely as possible with the quantities 
of the corresponding parent quark. This can be achieved with jet corrections as described in Jll|, |l2[]: after jets are 
reconstructed with the UA1 cone (Rc = 0.4) algorithm, their quality is improved in two steps. Firstly, single jets 
are corrected by taking into account out of cone energy and extra energy which does not come from the original 
parton. This energy correction improves the Higgs mass resolution by 35%. In the second step, the final state 
gluon radiation is corrected (FSR correction). Herefore two jets are combined into one jet taking into account 
geometrical effects. The second step gives an additional 10% improvement. Clearly, the single jet corrections 
are more important than FSR corrections, but nevertheless the FSR corrections can still improve the Higgs mass 
resolution significantly. Further improvements, not studied here, are possible by using tracks. 
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A detailed description of the analysis follows: 
o Trigger and Preselection 

Events are triggered, if there is one isolated lepton with px > 20 GeV/ c in the tracker acceptance. A lepton is 
considered as isolated, if there are no additional tracks with pr > 2 GeV /cm a cone of 0.3 around the lepton. The 
lepton reconstruction efficiency is assumed to be 90%. The preselection includes the requirement of at least two 
jets with E T > 30 GeV and \r)\ < 2.5. 

o Pile Up Jet Rejection 

In order to be able to count the jets of one particular event, pile up jets are removed from the event: jets are 
removed, if they contain no track with p? > 2 GeV jc in a cone of 0.4 around the jet axis, or if they contain a 
track of p? > 2 GeV jc in a cone of 0.2 around the jet axis which is coming from the wrong primary vertex: 
A(V^ rack , V%) < 250 fun. The identification of primary vertices is possible, because the isolated lepton points to 
the correct primary vertex and the primary vertex z coordinate (VS) resolution is as good as 25 /im, whereas the 
z bunch crossing spread has a = 53 mm, as illustrated in Figure || The performance of this procedure depends 
basically on instantaneous luminosity (number of pile up events) and V% resolution of the CMS tracker. 
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Figure 3: Example of bunch crossing spread in z direction: 21 spikes indicate the primary pp collision vertices due 
to one hard interaction (marked with a triangle) and twenty superimposed pile up events. The bin width of 250 p/m 
corresponds to ten times the resolution of the z coordinate of the primary vertex V% . The shaded area is a measure 
of the expected frequency distribution of primary vertices. 



o b-tagging 

Two tagged jets are required per event. A jet is considered as &-tagged, if there are two tracks with impact parameter 
significance cr(ip) > 2.0, or if there are three tracks with a(ip) > 1.6, or if there are two tracks with u(ip) > 1.0 
and a lepton inside the jet. Here ip is the signed transverse impact parameter of a track with pt > 0.9 GeV jc and 
ip < 0.2 cm. 

o Jet Corrections 

The tagged jets are corrected as fe-jets with or without lepton and additional jets which are closer than 0.6 to one of 
these jets are corrected as gluon jets and combined with the 6-jets according to the procedure described in [pTl n2|. 



o Jet and Lepton Veto 

If there is an additional jet with Et > 20 GeV and |^| < 2.5, the event is rejected. Also events with an additional 
isolated lepton (pr > 5 GeV jc) are not accepted. The jet veto is only possible after "pile up jet rejection". 

o Relative E T Balance 

The relative Et balance is defined as the ratio of the transverse energy of the event (Et of b, b, I and v which is 
ideally close to zero) and the Et sum of the four objects expected in the event. 

E ^ l M <0 .15 Er(u)°*EF" 



E T {b) + E T {b) + E T (l) + E T {y) 



O W^ 1 Identification The reconstruction of the transverse mass of the W boson is intended to suppress further 
badly reconstructed events. It turns out that only an upper cut improves the signal visibility. 

m T {l,v) < 100 GeV/c 2 
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O Higgs Mass Reconstruction The invariant mass of both fully corrected fo-jets is calculated. Finally events are 
counted which satisfy fh ± 1.4 a with m and a from Table |2|. 



TUftO = 
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= 91.28 


95.67 


100.1 


104.5 


108.5 


113.5 


117.8 


122.0 


126.5 


131.2 


a - 


= 10.76 


11.04 


11.19 


11.90 


11.67 


11.66 


12.28 


12.68 


13.22 


13.73 



Table 2: Generated m#o, reconstructed m and a in units of GeV/c 2 . rrii nv (j, j) of the pure signal is fitted with a 
Gaussian to obtain the masses and widths used in fh ± 1.4 a. 

The main background rejection comes from fe-tagging. The non-fo-jet backgrounds (mainly W^jj) are reduced 
strongly, but backgrounds with two genuine 6-jets are not suppressed. A large fraction of the top background can 
be removed with the combination of "jet and lepton veto", "relative Et balance" and "W identification". These 
steps also help to remove events with a large amount of unclustered energy and badly reconstructed events. All 
backgrounds are reduced by the mass window cut which is more effective with improved Higgs mass resolution. 




Figure 4: W^H signal (white or dark shaded) plus resonant W ± Z° background (shaded) plus non resonant 
background (light shaded). The crosses are signal plus background with statistical error bars and are well above 
the background in the signal region. Simulated Higgs mass is mfjo — 115 GeV/c 2 with Li nt = 300 fb~ x . 



For an integrated luminosity of 300 fb \ 1610 W^H events (m H o =115 GeV/c 2 ) are selected with an effi- 
ciency of 5%, starting from all triggered events. 1 198 W^Z events are selected with an efficiency of 0.3%. The 
efficiency for W ± jj events (27565) is 0.003%, for tt events (36089) is 0.1% and for tb events (6096) is 0.09%. 
Figure Q shows the expected signal plus background. The signal peak can be hardly seen. If one considers the 
simulated data points with their statistical errors, there is however a statistically significant signal above the back- 
ground. Background subtraction can improve the visibility and is done here in two steps. First the non resonant 
background is fitted with a polynomial of degree eight and then subtracted. It is important to know the shape of 
this background. This information can be obtained from detailed Monte Carlo simulations or experimentally by 
varying the 6-tagging quality. Figure |5] (left) shows the result of this subtraction: one can see a double peak which 
is a superposition of the Higgs boson and the Z° peaks. Even for a very small Higgs signal the Z° peak should be 
still visible. This is a good cross check of the subtraction method, and can be used to calibrate the Higgs mass. In a 
second step the W^Z events are subtracted. The accurate magnitude of this background can be estimated easily 
from the analysis of W^Z — > l^vl + l~ events. The result is shown in Figure |5] (right) and the Gaussian fit of the 
Higgs peak is in good agreement with the simulated pure signal and the Higgs mass can be fitted. 
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Figure 5: Left: signal plus background after subtraction of non resonant background. The simulated data points 
are fitted with two Gaussians. Right: signal plus background after subtraction of all backgrounds. The simulated 
data points are fitted with a Gaussian which describes the pure W ± H° signal (shaded) well. 



With this type of analysis we obtain following results for mjjo =115 GeV/c 2 and Li n t = 300 fb lm . signal to 



background ratio S/B = 2.3% , significance S/yB = 6.0 , precision on WWH coupling Agwwn /gwwn = 
8.4% and precision on the mass Am/m — 2.3%. S and B are the number of events in the mass window around 
the Higgs mass peak. The precision on the mass is obtained from the Gaussian fit in Figure pi (right). 



3 Expectations for SM Higgs 

The SM results are shown in Figures ^ and 0. For L int = 300 fb" 1 the 5 a discovery limit is at TOjyo < 123 
GeV/c 2 . A discovery during the low luminosity phase (L int < 30 fb" 1 ) is not expected in this channel - even in 
the easiest cases i.e. at low Higgs masses already excluded by LEP. Higgs boson discovery is expected in another 
channel. The W ± H° channel is nonetheless interesting, because it allows the measurement of the Higgs couplings 
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Figure 6: S/B, S/ \/B (Lint — 300 fb 1 ) and Li n t required for S/vB = 5 versus generated SM Higgs mass in 
the W^H — + V^vbb channel. No k-factors (k= 1.0) are used. 
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which will tell us something about the nature of this particle. The WWH coupling is proportional to the square 
root of the number of signal events. The expected precision in the measurement of the WWH coupling is given 
in Figure [7| and is of the order of 10% assuming a known branching ratio BR(H° — > 66). 
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Figure 7: Precision of WWH coupling as a function of generated Higgs mass in the SM. Li nt = 300 fb 1 is 
assumed and no k-factors are applied. 



Even if it is impossible to determine BR(H° — ► bb) in absolute value, the different signals of various Higgs 
production mechanisms can be exploited to determine ratios of Higgs couplings. For example, the ratio of gwwH 
and y t (top Higgs Yukawa coupling constant) is independent of BR(H° — > bb) and can be determined with the 
tiH° and W ± H° channels. The ratio of BR(H° -► 66) and BR(H° -> 77) can be determined by comparing 
W ± H°, H° — ► 66 with l^if , 7?° — * 77 final states produced by the same qqi — * iy* — > W ± H° mechanism. 



4 MSSM Results 

To give an idea about the discovery potential of the corresponding channel W^hP — > V^vbb in the MSSM, 
we extrapolate the SM results (by rescaling the production cross section times branching ratio, obtained with 
SPYTHIA flrjj]) and discuss the parameter space coverage of one benchmark scenario called maximum m-h sce- 
nario [53]. This scenario turns out to be the most unfavourable one as h° — > 66 visibility decreases with increasing 
/i° mass. The reason is the rapidly falling cross section and branching ratio with increasing Higgs mass which 
limits the discovery potential of this channel in the SM as well. 
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Figure 8: Discovery contours in the MSSM (maximum mh scenario) parameter space for L{ n t = 300 /6 _1 (left) 
and for Li nt — 600 fb^ 1 (right). S/y/~B > 5 to the shaded side of the solid line. The dotted and dashed lines are 
the isomass curves form h o = 125 GeV/c 2 andm h o = 115 GeV/c 2 , respectively. 



6 



Figure |8| shows the parameter space coverage in the itia tan [3 plane for two integrated luminosities integrated 
over very long LHC running periods (5 - 10 years). In both cases there is an inaccessible region at low itia', the 
second difficult region at high rriA and tan (3 disappears with increased integrated luminosity. In scenarios with 
non maximum mu the difficult regions are somewhat reduced implying that for high enough integrated luminosity 
most of the MSSM parameter space can be explored with this channel. 



5 Conclusions 

At this level of simulation, including high luminosity event pile up, we conclude that the selection and reconstruc- 
tion of the W^H — > l^vbb signal is possible and provides useful information. Excellent 6-tagging performance 
and good mass resolution is crucial for a successful analysis. In addition, track and primary vertex reconstruction 
is important for the separation of jets from the hard interaction from jets of pile up events. 

A very high integrated luminosity of 300 fb^ 1 is necessary in order to have reasonable sensitivity to this signal: in 
the SM the discovery is limited at a Higgs mass of < 123 GeV/c 2 . In the MSSM a large fraction of the parameter 
space can be covered. It is clear that the very large integrated luminosity required makes this channel inappropriate 
for a first discovery of the Higgs boson, in contrast to the channel UH°, H° — > bb jl5[]. 

The W^H — > l^vbb channel allows direct measurement of the WWH coupling for known branching ratio 
BR(H° -> bb). Even if BR(H° -> bb) is not known, this channel can be combined with other channels with 
H° — > bb decay [^|] and ratios of the corresponding coupling strengths can be determined which would help to 
better understand the nature of the Higgs boson. 
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